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Abstract 

We test the hypothesis of non-radiative neutrinos decay using the latest IceCube data. Namely, we calculate the 
track-to-shower ratio expected in IceCube for the normal and inverted neutrino mass hierarchy taking into account the 
uncertainties in neutrino oscillation parameters. We show that the subset of data with energy above 60 TeV actually 
excludes the possibility of a neutrinos decay at the 1 sigma level of significance for both neutrino mass hierarchies. 
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1. Introduction 

The neutrinos decay scenario has been early proposed 
as an explanation of the solar neutrinos problemlT] 13. 
Although the neutrinos decay through radiative pro¬ 
cesses is well constrained a, the lifetime limits for 
non-radiative channels are very weak. Generally re¬ 
searchers focused their attention to non-radiative pro¬ 
cesses of the kind v' —> v + X, where neutrinos de¬ 
cay into possibly detectable neutrinos (or antineutrinos) 
plus truly invisible particles, X, e.g., light scalar or pseu¬ 
doscalar bosons. As an example, the neutrino decay can 
take place through the coupling of the neutrino to a very 
light or massless particle, such as a Majoron, which is 
also responsible for spontaneous neutrino mass genera¬ 
tion ||ll|5l. 

In general, considering a neutrinos source, decay 
could deplete the flux of neutrinos by the factor. 
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where m, and r, are the masses and lifetime of the i-th 
neutrino, subject to decay, E is their energy and T is the 
time since production. As a consequence, the sensitiv¬ 
ity to the unknown ratio t/ot changes for different dis¬ 
tances and energies. Actually the strongest reliable limit 
for non-radiative processes and for hierarchical masses, 
namely r/m > 1.1 ■ 10“^ s/eV, is obtained by the non¬ 
observation of solar electron antineutrinos in Kamland 
||6l . This limit can be in principle improved by sev¬ 
eral order of magnitude by using the high energy cosmic 
neutrinos observed in detectors such as IceCube Q. 

This intriguing possibility became more interesting 
on the light of the recent search for High Energy Start¬ 
ing Events (HESE) in IceCube detector that provided 
the first evidence for high-energy cosmic neutrinos|[8l 
|9][T3. In three year of data taking ISl, 37 events with 
deposited energies above 30 TeV were observed, rela¬ 
tive to an expected background of 8.4 + 4.2 cosmic ray 
muon events and 6.6+5.9 atmospheric neutrinos. 

Recently, Palladino et al. lfTTI have discussed the 
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compatibility of these data with the hypothesis of cos¬ 
mic origin. In order to reduce the background contam¬ 
ination, authors have considered the subset of events 
with deposited energy above 60 TeV. Using the addi¬ 
tional information provided by muon neutrinos passing 
through the Earth, they show that the observed track-to- 
shower ratio matches expectations for neutrinos of cos¬ 
mic origin. 

In this Letter, we discuss the impact of the non- 
radiative neutrino decay scenario on the expected Ice- 
Cube signal in term of the track-to-shower ratio. More 
precisely, considering the actual uncertainties on oscil¬ 
lation parameters, mass hierarchy and production mech¬ 
anisms; we use the observed events to constrain the hy¬ 
pothesis of neutrinos decay. As a remarkable result, we 
show that the subset of data with energy above 60 TeV 
excludes the possibility of a neutrinos decay at the 1 
sigma level of significance for both neutrino mass hier¬ 
archies. 


2. Complete Decay Scenario 


To maximize the effect of neutrino decay, we consider 
the phenomelogy of a generic neutrinos decay and we 
assume that decays are always complete (i.e., that the 
exponential factors vanish). This last assumption is rea¬ 
sonable for Cosmic Sources considering that L ~ \ Gpc 
and £ ~ 100 TeV), namely 


TjE 


10 ^ 


( 1 ) 


where we used L - cT. We call this hypothesis Com¬ 
plete Decay Scenario. In this case, the flux of neutrinos 
of flavor { expected to the Earth is; 


t' j=stable 


where U is the neutrino mixing matrix and the sum only 
involves stable mass eigenstates. Eollowing iTTi nTII . we 
introduce the flavor fractions at Earth (i.e., in the detec¬ 
tion point), defined as: 

= O^/chtot with (htot = Y 

£=acTive 

In term of this quantity the relation between the flavor 
ratios at Earth and the ones at source becomes: 

where the indeces j and E run over stable mass eigen¬ 
states and active flavors, respectively. The expression 


in Eq.([^ also holds in the case where sterile neutrino 
states are included. Different production mechanisms 
generate a different initial flavor content at the source 

^ : ^). Eor example the most studied production 
processes are; 

(1/3 : 2/3 : 0)for n decay; 

(1/2 : 1/2 : 0) for charmed mesons decay; 

(1 ; 0 ; 0) foryS decay of neutrons', 

(0:1:0) for tt decay with damped muons. 

Eor cosmic neutrinos, without neutrino decay processes, 
the final flavor ratios are expected to be very near to 
: ^t) ~ (1/3 : 1/3 : 1/3), independently on the 
specific production mechanism KElinilllinillll. On 
the other hand, turning on the ’’invisible” decay process 
discussed above, the flavor ratios at the Earth change 
dramatically Q. In the case of only one stable mass 
eigenstate j and for the standard scenario with three ac¬ 
tives neutrinos, the expression of Eq.([^ becomes 

& = \Uij\\ (4) 

where the sensitivity to the different production mech¬ 
anisms has disappeared. The simplified expression in 
Eq.Q is due to the fact that only one mass eigenstate is 
stable and all the others have had time to decay, so that 
the universe is populated only by neutrinos in the state 
of mass j and obviously the content of flavor { expected 
to the detector is \ U(j\^. 

Assuming normal mass hierarchy, the simplest case 
is to consider that both V3 and V2 decay, so that only 
the lightest eigenstate vi survives. In this case the fi¬ 
nal flavor ratios are expected to be (|t/<.i|^ : |t/yjip : 
Wrif') ~ (0.68 : 0.21 ; 0.11). Eor inverted mass hi¬ 
erarchy, a similar expression holds but the stable state 
in this case is V3, i.e. the final flavor ratios are expected 
to be (\Ue 3 \^ : \U^ 3 \^ : |t/^ 3 p) ~ (0.02 : 0.57 : 0.41). In 
order to evaluate these ratios, we used the best-fit values 
for the oscillation parameters reported in IfTTl . 

3. Track-to-Shower Ratio 

To compare the expectations with the data collected 
by IceCube we need to convert the flavor neutrino ra¬ 
tio expected to the Earth in the corresponding observ¬ 
able quantity, i.e. the track-to-shower ratio including 
the detector response. Indeed, by exploiting the event 
topology, the flavor discrimination is possible, in prin¬ 
ciple, in IceCube detector. The HESE data collected 
by IceCube during T = 988 days, encompass two dif¬ 
ferent topologies: ‘shower’ topology, that includes neu¬ 
tral current (NC) interactions of all neutrino flavors and 
charged current (CC) interactions of and vp, ‘track’ 
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Figure 1: Expected track-to-shower ratio for decaying cosmic neutrinos. The distributions show the effect of uncertainties in the 
neutrino oscillation parameters. The left (resp. right) distribution is obtained for normal (resp. inverse) hierarchy. The black line 
shows the normalized likelihood function of the IceCube data-set discussed in the text. 


topology produced by CC interactions of v^. Thus the 
track-to-shower ratio is the observable quantity that is 
most directly related to the flavor ratio at the Earth as 
discussed in several recent papers iniiisi [211201 EH. 

In particular Palladino et al. lITTI described the way to 
calculate the number of showers and tracks N-y in the 
IceCube detector for an isotropic flux Of of neutrinos 
and antineutrinos described as 


effective area A^(£’) giving tracks, this fraction is about 
~ 0.8 and mildly dependent on energy as discussed in 
na. So that, the following useful expression to obtain 
the track-to-shower ratio from the final neutrino flavor 
ratios can be obtained generally: 

^ -, ( 8 ) 

Ns C\ C2^^i + Cs 


Of(£) = 


Fe ■ 10 -^ 
cm^ s sr GeV 



(5) 


where the factors F( are (non-negative) adimensional 
coefficients and a is the spectral index. The number of 
showers and tracks can be obtained by using the follow¬ 
ing equations: 

Ns ^ T \^Fe Ce -i-FrCr+ F C^(l - pj)} , (6) 

Nt = TFi^c^pt, ( 7 ) 


where T is the time of exposure and the coefficients ci - 
At: J dE Af(E) include the detector effective areas 
Af(£’)lT0l and the energy dependence of the fluxes. The 
extra factor px in Eq.(0 accounts for the fraction of the 


in which cj are given, with good approximation, by: 

Cl = ^ 2.23 -H 1.15(a - 2) + 0.46(a - if 

C2 = == -1.99 - 1.15(0' - 2) - 0.46(a - 2)^ 

C 3 = ^ -0.56 - 0.78(0' - 2) - 0.40(o' - if 

In the context of neutrino decay scenario described 
above and fixing the spectral index a = 2, we have that 
for Normal mass hierarchy Eq. becomes: 

^ ^ _ .g. 

Ns 2.3-2.0|t/j2-0.6|t/,/’ 

where the index j indicates the stable mass eigenstate. 
Eor Normal Mass Hierarchy j - I and substituting the 
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best-fit values of the oscillation parameters we obtain 
= 0.12. For Inverted Mass Hierarchy y = 3 and 
this fraction becomes - 0.62. The impact of the 

spectral index of neutrino spectrum is not so dramatic, 
in fact if we consider a greater slope, like a - 2 . 6 , these 
number become = 0.09 and - 0.56. In¬ 

stead, a much more important contribution is given by 
uncertainties on neutrino oscillations parameters, that 
allow these numbers to fluctuate in a more large interval. 
To correctly account for these uncertainties we construct 
likelihood distributions of sin^ 612 , sin^ 613 , sin^ 023 and 
6 from the profiles given by ini. Namely, we as¬ 
sume that the probability distributions of each parame¬ 
ter are provided by X = exp (-A;^f^/2). Then, we com¬ 
bine the various likelihood functions assuming negligi¬ 
ble correlations and we determine the probability distri¬ 
butions of Nt/Ns by MonteCarlo extraction of the os¬ 
cillation parameters. The resulting Probability Density 
(PD) distributions are reported in Fig.Q with a color 
code blue (resp. red) for Normal (resp. Inverted) Mass 
Hierachy. 

The predicted PDs can be compared with the distri¬ 
bution of the track-to-shower ratio preferred by the Ice- 
Cube data. In the first column of Tab. ([T]) we report the 
events observed by IceCube in the HESE data set ||8l 
with a deposited energy above 60 TeV. We devided the 
events by topology and we includ, in brackets, the re¬ 
spective number of estimated background events due to 
atmospheric neutrinos and muons. 

The number of tracks Nj and showers Ns which have 
to be ascribed to extraterrestrial sources can be esti¬ 
mated from the Poisson likelihood functions: 

X(M) ex df X g-'*' (10) 

where di = W H- and the index i = T, S is used to 
refer to track and shower events and are the observed 
events. We are assuming that the prompt atmospheric 
neutrinos give negligible contributions, as it required by 
the spectral and arrival angles distributions of IceCube 
events. A completely equivalent and independent infor¬ 
mation can be obtained by the recently released IceCube 
data on passing muons 1221 . as already done by Pal- 
ladino et al. in na. In the assumption of £ ^ neutrino 


Table 1: The total number of events, ni, with the expected background, 
bj, seen at IceCube, in the energy range 60TeV < E^ep < 3PeV. 



Hi {bi) 988 days 

(bj) 1460 days 

Track 

4(2.1) 

8(3.1) 

Shower 

16 (0.7) 

25 (1.0) 

Total 

20 ( 2 . 8 ) 

33 (4.1) 


spectrum for the passing muons, the flux normalization 
corresponds to = 1.01 +0.35 and taking into account 
the equivalence between £), and Nj expressed by the 
following equation (in which T is the exposure time): 

Ns = (2.94 xFe + 0.33 x + 2.20 x F^) ■ T 
Nt = 1.31x£^-7’ 

we can include also this information in our analysis by 
constructing a combined likelihood, given by the prod¬ 
uct of the 2 Poisson likelihoods for Nj and Ns and of 
the Gaussian likelihood for F^ (i.e. Nj)- Then we ex¬ 
tract the bounds on the track-to.shower ratios of events 
ascribed to cosmic neutrinos by marginalizing with re¬ 
spect to the total number of events. We obtain: 

[988 days] ( 11 ) 

Ns 

where the error was obtained by integrating out sym¬ 
metrically (1 - CL)/2 on both sides of the Nt/Ns distri¬ 
bution using a confidence level CL = 68.3%. 

The normalized likelihood distribution of the HESE 
data is reported in Fig. 0 with a black line and the 1 
sigma region is also indicated to easily see its sovrap- 
position with the predicted PDs in the neutrino decay 
scenario. 

4. Discussion and conclusion 

Repeating the analysis with 4 years data, under the 
hypothesis that the background is simply proportional 
to the exposure time, we found a new range given by: 

^=0.20ti° [4 years] (12) 

We can notice that this interval is compatible with that 
found for the 3 vear liTTI also if the uncertainties are 
slightly smaller with new data; the track-to-shower of 
4 years is more favorable for pions decay at the source, 
respect to that obtained with 3 years data. In conclusion, 
it is important to remark that the non-radiative neutrino 
decays are excluded at least at 1 sigma for both neutrino 
hierarchies and for both sets of observed data. Moreover 
these are general considerations and have a poor depen¬ 
dence from the slope of the spectrum of neutrinos. 
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